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The renormalization group technique is applied to one-dimensional electron-phonon Hubbard 
models at half-hlling and zero temperature. For the Holstein-Hubbard model, the results of one-loop 
calculations are congruent with the phase diagram obtained by quantum Monte Carlo simulations 
in the {U,gph) plane for the phonon-mediated interaction ^ph and the Coulomb interaction U. The 
incursion of an intermediate phase between a fully gapped charge-density-wave state and a Mott 
antiferromagnet is supported along with the growth of its size with the molecular phonon frequency 
(jUq. We hnd additional phases enfolding the base boundary of the intermediate phase. A Luttinger 
liquid line is found below some critical f/* followed at larger U ~ ^ph by a narrow region of 

bond-order-wave ordering which is either charge or spin gapped depending on U. For the Peierls- 
Hubbard model, the region of the {U,gph) plane with a fully gapped Peierls-bond-order-wave state 
shows a growing domination over the Mott gapped antiferromagnet as the Debye frequency ud 
decreases. A power law dependence gph ~ is found to map out the boundary between the two 
phases, whose exponent is in good agreement with the existing quantum Monte Carlo simulations 
performed when a hnite nearest-neighbor repulsion term V is added to the Hubbard interaction. 

PACS numbers: 71.10.Fd, 71.30.+h, 71.45.Lr 


I. INTRODUCTION 

In approaching the physics of highly correlated electron 
systems, we often come against experimental evidence 
for significant coupling between interacting electrons and 
lattice degrees of freedom of different kindJ^lt^. The fact 
that the electron-phonon coupling can play an impor¬ 
tant part besides the electron-electron interaction takes 
on particular importance in reduced dimensions where 
the intrinsically expanded range of electronic instabili¬ 
ties can be further broadened by the interplay between 
both interactions.®® 

In this matter one-dimensional (ID) models of inter¬ 
acting electrons coupled to a bosonic phonon field have 
been long considered as primary models to study the 
competing effects of Coulomb and retarded interactions 
on ordered phases at zero temperature. In the present 
work we shall be concerned with two models at half- 
filling that have been the focus of considerable atten¬ 
tion in the past few decades, namely the ID Holstein- 
Hubbard model (HHjzE 

and the Su-Schrieffer-Heeger- 
Hubbard modeP and its variant for optical phono ns 
called for short the Peierls-Hubbard (PH) modeP®^. 

The corresponding Hamiltonians can be introduced 
through the following common form in Fourier space, 

H=Y, s{k)cl^Cka + + \) 

k,a q 

+ X 3ik, q){bl + 

^ k,q,(j 

+ ^ ^ Vq ( 1 ) 

q q 


where c|.^(c/ccr) creates (annihilates) an electron of wave 
vector k and spin a, and L is the number of lattice sites 
(the lattice distance a = 1). The electron spectrum is 

e{k) = —2tcosk, (2) 

where t is the hopping amplitude. At half-filling the band 
is filled up to the Fermi points ±kF = ±7r/2. 

In the case of the HH model, momentum independent 
intramolecular phonons of energy ujq = ujq {h = 1), de¬ 
scribed by the creation (annihilation) phonon operators 
6j(6g) of wave vector g, are coupled to electrons through 
the momentum independent coupling constant 

g{k,q)^go^ (HH) (3) 

where > 0- 

For the Su-Schrieffer-Heeger model, the phonons be¬ 
long to an acoustic branch of frequency, 

LUq = cjdI sing/2|, (4) 

where ujd corresponds to the Debye frequency for the 
zone boundary q = ir phonon, namely at twice the Fermi 
wave vector 2k f for a half-filled band. Their coupling 
to electrons results from the modulation of the electron 
transfer integral by lattice vibrations which leads to the 
momentum-dependent coupling constant 

5(^,9) = *5Dsin|cos (A:+I) (5) 

where is a constant proportional to the spatial deriva¬ 
tive of the hopping integral t. Retaining only the coupling 
to q = 2kF phonons at the zone edge corresponds to the 
PH limit. 
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The Coulomb part of the Hamiltonian Q is common 
to both models. It comprises the Hubbard and extended- 
Hubbard interactions U and V {Vq = V cos q) between 
electrons on each site and nearest-neighbor sites. These 
are expressed in terms of the occupation number, Uq = 

ng-f + Ugi, where Ug^ = A+g^a^^ka- 

The phase diagrams of both models have been the ob¬ 
ject of sustained interest during the last decades. To start 
with the half-filled HH model in the pure Hubbard limit 
at H = 0 , many efforts have been devoted to clarifying 
its structure. Special emphasis has been put on the fron¬ 
tier separating the fully gapped (spin and charge) 2k f 
site-centered charge-density-wave (CDW) and the charge 
(Mott) gapped 2kp spin-density-wave (M-SDW) phases, 
namely when the amplitude of the phonon-mediated elec¬ 
tron coupling 

^ph — 

becomes of the order of the local repulsion U. 

The existence of an intermediate metallic state was 
predicted to occur long ago on the basis of real-space 
renormalization group argumentJ^^, but found to be ab¬ 
sent in former world line Monte Carlo simu lation s^ and 
two-cut off renormalization group analysiS^^^. The 
presence of such a state was more recently supported 
by density matrix renor malization group (DMRG) and 
variational methodS^^^. Its origin was thereafter exten¬ 
sively discussed from various approache^^^. In par¬ 
ticular, quantum Monte Carlo (QMC)^* ^^ * ^^ * and fur¬ 
ther dmrgP calculations have corroborated its exis¬ 
tence over a definite, cjo—dependent region of the (I/, ^ph) 
plane. On more analytical grounds, the phase diagram of 
the HH model has been investigated W the functional RG 
(fRG) method at the one-loop level.^ The frequency de¬ 
pendence of the electron-electron vertices has thus been 
obtained in the putative intermediate region, umklapp 
scattering, though irrelevant in the static limit, was found 
to be large at finite frequency. This was interpreted as the 
driving force of a GDW phase gapped in the charge sec¬ 
tor, going against the existence of an intermediate metal¬ 
lic state, or at the very least the existence of dominant 
superconducting correlations whose presence in numer i- 
cal simulations was attributed to finite-size effect^^^^. 

In the first part of this paper, we reexamine this prob¬ 
lem using a RG method^ similar to the fRG and pro¬ 
ceed to a detailed scan of the phase diagram of the HH 
model as extracted from the singularities of susceptibil¬ 
ities and the frequency dependence of electron-electron 
vertices. We thus reproduce to a high degree of accu¬ 
racy the boundaries of the intermediate phase found by 
the QMG simulations of Hardikar et al^ at ^fph ^ U 
and arbitrary cJcP^. In this sector of the phase diagram, 
the GDW and to a lesser degree singlet superconducting 
(SS) response functions are found to be singular. This is 
compatible with a metallic state characterized by effec¬ 
tive attractive couplings, which evolves toward those of 
the attractive Hubbard model at half-filling when the non 
adiabatic, cjo ^ oo limit is taken at 1/ = 0. The present 


RG calculations further predict the presence of a particu¬ 
lar structure of the boundary between the metallic phase 
and the M-SDW state. A gapless Luttinger liquid and a 
gapped 2kF bond-centered charge-density-wave (BOW) 
phases are successively found to enfold the boundary as 
a function of U. Both prevent a direct passage from M- 
SDW to GDW in the (I/,^ph)—plane at finite ujq. 

As to the phase diagram of the PH model at zero tem¬ 
perature, it is known to display a simpler structure. Nu¬ 
merical simulations of the phase diagram from QMCP^ 
or DMRGP^ techniques for instance agree with a direct 
transition line between the fully gapped Peierls BOW (P- 
BOW) state and the M-SDW phase in the (I/, ^fph)—plane 
for small R > 0, a feature of PH phase diagram also 
well depicted a nalytica lly by the the so-called two-cutoff 
RG calculation d^^ l ^^ l ^ , provided the couplings and De¬ 
bye frequency ujd are not too large. The RG calculations 
to be developed here in the continuum limit agree with 
QMG results carried out for non dispersive phonons in 
the weak coupling sector for U and 

In Sec. H, one-loop flow equations for the scattering 
amplitudes, self-energy and response functions are given 
for both models in the electron gas continuum limit. In 
Sec. HI, the flow equations are solved for the HH and PH 
models, respectively. The singularities in the couplings 
at zero and finite frequency, along with static suscepti¬ 
bilities are tracked in the (I/,^ph)—plane and serve to 
map out the phase diagrams for different and V in 
the case of the PH model. Gomparison with numerical 
results is made. We conclude in Sec. IV. 


II. RENORMALIZATION GROUP EQUATIONS 
OF THE CONTINUUM LIMIT 

From the functional integral formulation of the parti¬ 
tion function, the integration of the phonon field leads to 
a (Matsubara) frequency dependent interaction between 
electronJ^. When combined with the Goulomb terms U 
and V, this yields an effective electron-electron interac¬ 
tion of the form: 


^i=l,2,3(<^l5 ^2, ^ 3 ) 


ffph,i 

1 + (oji F 


(6) 


when expressed in the ^-ology picture of the electron gas 
model, where uji = (2ni + l) 7 rT. In this continuum frame¬ 
work, the electron spectrum 0 is linearized around the 
Fermi points pkp = -^kp and becomes 


e{k) ^ ep{k) = vpipk — kp)- (7) 

The spectrum is bounded by the band width cutoff, 
EqI 2 = pvpl2^ on either side of the Fermi level, where 
vp = 2t the Fermi velocity at half-filling. The interac¬ 
tion defined on the Fermi points then breaks into three 
pieces, namely the backward (^ 1 ), forward (^ 2 ) and umk¬ 
lapp (^ 3 ) scattering amplitudes whose frequencies satisfy 
cji -f u;2 = CJ3 -h CJ4. 
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The corresponding attractive amplitudes for the 
phonon-mediated part are given by 

^ph,l,2,3 = —^ph, (HH) (8) 

for the HH model, where ^ph = 2 ^q/ci;o- 

For the PH model, however, the electron phonon- 
matrix element is momentum dependent, which yields 

5^ph,l ^ph,3 ^ph5 ^ph,2 O 5 (^^)' (^) 

where 

9ph = 2g%lujD- (PH) 

Therefore repulsive umklapp and the absence of forward 
scattering at the bare level are direct consequences of 
the momentum dependence of the bond electron-phonon 
interaction <§• This is a key ingredient for the differences 
between HH and PH models. 

For the Coulomb part that is common to both models, 
we have the familiar bar e coup ling combinations of the 
electron gas at half-filling 

gi=g 3 = U-2V, g 2 = U + 2V. (10) 


All the bare amplitudes are confined to the weak-coupling 
regions Qi = gij'KVF < 1 and ^ph = < 1- 


We apply a Kadan off-Wilson RG procedure for the par¬ 
tition fnnctionj ^^ l ^^ l which consists at zero temperature 
of the successive integration of electronic momentum de¬ 
grees of freedom from the cutoff energy ±F^o/2 down to 
the energy distance ±£’o('^)/2 on either side of the Fermi 
level. Here — Eoe~^ is the effective bandwidth at 

step £. This successive integration in the momentum is 
performed for all Matsubara frequencies. The frequency 
axis is divided into 61 sections or patches between the 
cutoff values TcJmax = ±1.5£’o/2, which serve as bounds 
of integration for the frequency. The interactions are 
taken constant over each patch where the loop integrals 
are done exactly. The successive integration leads to the 
flow of coupling constants, one-particle self-energy, and 
susceptibilities. At the one-loop level, the flow of the 
normalized scattering amplitudes at T ^ 0 reads 


degi — 2 


9^^2 = 


(cji,u;2,cj3) = 2 


/ I [ - gi{uji,uj\ujs)gi{uj4,uj',uj2) - ^ 3 (^ 1 , ^3)^3(^4, ^ 2 ) +^i(cJi,Cj',CJ3)^2(^',^4,^2) 

J uj' ^ 

53 (^ 1 , w',W3)53(w',W4,W2)]/+(w',Wi - W 3 ) - 51 (wi, W 2 , w')g 2 (^3, W 4 , w')/-(w', +W 2 )|, (H) 

j I [ - 5 l (‘^l, ^^ 2 , w') 5 i(w 3 , W4, w') - 52(^1, W 2 , w') 52 (w 3 , W4, w')] J_ {u}', Wi + W2) 

+ [52(w',W2,W3)g2(w',W4,Wi) + 53 (w', W 2 , W 3 )g 3 (w', W 4 , Wi)]/+(w', W 2 - W 3 )|, ( 12 ) 

/ I [ - 25i(wi,w',W3)g3(w4,w',W2) +5i(wi,w',W3)g3(w',W4,W2) 

J Uj' ^ 

+ 52 (w',W 4 ,W 2 ) 53 (wi,a;',W 3 )]J+(w',a;i - 023) + 52(02', W2,0)3)53(0;',024, wi)J+(w',022 - o; 3 )|j( 13 ) 


where di = djdi and duj' I27r. Here /=p refer 

to the 2kF electron-hole (Peierls) and electron-electron 
(Cooper) loops, 

/zp(cj, = 

{uj — X](cij))(ciJ ^ zb X](fl zb Ct’)) “b A^ 

[{uj - S(u;))2 + A2] [(CJ T ^ ± s(fl ± u;))2 + Kj] ’ 

(14) 

where A^ = Eo{£)/2. The imaginary part of the Matsub¬ 
ara self-energy, H, for right or left-going electrons obeys 
the equation 

deT,{Lo) = J | 5 i(o;',o;,o;i) - 252 ( 0 ;,o;',o;)| 


o;' - E(o;') 

"" "[o;'-S(o; 0]2 + Af 

with the initial condition i;(co’)|^=o = 0- The one-loop 
flow is essentially governed by the superimposition of 
Cooper and Peierls pairing channels which interfere with 
one another, mixing ladder, closed loop, and vertex cor¬ 
rection diagrams at every order. The combination con¬ 
trols the nature of correlation at long distance, and in 
the presence of a phonon part will be dependent on the 
degree of retardation or the phonon frequency u;o,d- 
The nature of correlations as a function of £ are an¬ 
alyzed from a selected set of static normalized suscepti¬ 
bilities Xfi{= TT'T’izx^), namely those likely to be singular 
in the density-wave and superconducting channels. The 
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static susceptibilities obey an equation of the form 

diXn = [ (16) 

J UJ 

Each of these involves a pair vertex function governed 
by the equation 

dez^{uj)= (17) 

J Uj' 

where z^{uj)\i=o = 1. For the charge density-wave chan¬ 
nel, we shall consider the CDW and BOW susceptibili¬ 
ties at the wave vector 2/ci?, and to which correspond the 
combinations of couplings 


g^{Lo\uj,Lo) 


= 52(w',w,a;) - 2 gi{uj,Lo\uj) 

^=CDW(BOW) 


± [53(w,w',w')-2g3(w',w,w')]. 


(18) 


For the spin-density-wave channel, the antiferromagnetic 
or site centered SDW susceptibility may be also singular; 
it corresponds to the combination of couplings 

5m=sdw(w',w,w) = g2{uj',oj,oj) (19) 

- [53(w,w',w')-253(w',w,w')]. 



£ £ 



£ £ 


In the superconducting channel, only the SS susceptibil¬ 
ity may develop a singularity. It is linked to the combi¬ 
nation, 

5ss(w',a;,w) = -§2(1^, - gi{-ui,u,u'). ( 20 ) 


The solution of equations as a function of £ for Xn 
the scattering amplitudes reveals the singularities in 
these allow the identification of dominant and subdomi¬ 
nant singular correlation^^. A divergence at a finite £c 
indicates the breakdown of the weak-coupling one-loop 
RG procedure. Nevertheless, it singles out an energy 
scale or a gap A = EFe~^^ at zero temperature that cor¬ 
relates with a gap in the spin or/and charge degrees of 
freedom at £c. The gap in the spin sector is associated 
with a singularity in the static, cj ^ 0, limit for the at¬ 
tractive backward-scattering amplitude (11), whereas in 
the charge sector, a gap is found for singular either static 
repulsive or attractive umklapp (13). 


III. RESULTS AND DISCUSSION 


A. Holstein-Hubbard model 

Let us first consider the HH model at U > 0. In 
Fig. [^typical £ dependencies for the susceptibilities are 
shown for a fixed bare phonon-mediated coupling am¬ 
plitude, ^fph = 2t, and phonon frequency ujq = t. At 
U = 0.5t, that is in the region well above the bisecting 
line ^ph = U, a divergence is found only for the CDW 


FIG. 1: (Color online) In (a)-(c), the flow of normalized 
susceptibilities Xi^ £ for the Holstein-Hubbard model for 
Qph/t = 2 at different values of U/t. In (d), the flow of suscep¬ 
tibilities on the boundary of spin-density-wave phase with the 
intermediate phase; insert: susceptibilities on the boundary 
with the charge-density-wave phase. In all cases the phonon 
frequency is fixed at Ct;o = L 


susceptibility [Fig. 
tracted from the 


ig. [T| (a 
Fi^( 


(a)]. The value of £c — 3.25 ex- 
J (a) leads to a one-loop CDW gap 
A = EFe~^^ that is sizably reduced compared to the 
adiabatic -mean field- valu^® Aq = of the 

pure phonon model at cjo ^ 0 and for which, £q = 7r/2 
at the coupling considered here. The value of A is also 
smaller but closer to the non adiabatic cjo ^ oc limit 
of Eqs. (ll[^) at U = 0, corresponding to an effective 
attractive Hubbard model, where Aoo = and 

^oo = TT < at the one-loop leveP^. Repulsive U, which 
primarily reduces double occupancy of electrons on sites 
decreases CDW correlations and in turn their coupling 
to 2kF molecular phonons through vertex corrections in 
Eq s. (pT ])-(p!^, can explain these downward shifts of the 
ga] 


Regarding now the structure of scattering amplitudes 
at finite frequency in this CDW phase, we follow Ref.l^^ 
and show in Fig. ^ (a) the contour plot of the projected 
^i=i, 2 , 3 (<^i 5 ^ 2 ,^ 3 j ia the (cji,ci; 2 )-plane when £ ^ £c. As 
previously noted in the framework of fRCt^, the scatter¬ 
ing amplitudes display significant structure as a function 
of frequency, notably an attractive singular behavior for 
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FIG. 2: (Color online) Contour plot of the frequency depen¬ 
dent normalized couplings gi(cui in a finite portion of 
the (ct;i,ct; 2 ) plane at ^ph = 2t, cjq = t, ^ close to the critical 
ic, and for two values of Coulomb interaction: (a) U = 0.5t 
(CDW) and (b) U = At (M-SDW). 


the CDW phase that persists away from the origin of the 
(cji,ci; 2 )-plane. Divergences toward negative values are 
seen in both CJ 2 ,cji) and with max¬ 

ima at the origin [Fig.|^(a)], supporting the existence of 
a gap in both spin and charge degrees of freedom compat¬ 
ible with the absence of enhancement in spin correlations 
and the existence of commensurate CDW order at half¬ 
filling. 

One can wonder about the impact of vertex singulari¬ 
ties in the finite frequency range on correlations. In the 
present RG procedure both the marginal and irrelevant 
(frequency) parts of couplings are comprised in the three- 
variables scattering amplitudes ^'^(co’i,co’ 2 ,cjs). However, 
one can get a rough idea of how the frequency depen¬ 
dence emerges. Assuming that the amplitudes can be 
expanded around the origin as 


3 

W2, W3) = gj(0, 0, 0) + + . . . , 

i=i 

we see that the amplitudes, gij, etc., of the frequency 
dependent part are by power counting superficially ir¬ 
relevant in the RG sense. However, these will be cou¬ 
pled to the marginally relevant ^^(0,0,0) at the one- 
loop level and may become in turn relevant when one 


or more ^^(0,0,0) scales to strong coupling, that is as i 
approaches £c- This is likely what happens in the present 
RG procedure. Being strongly weakened by the fre¬ 
quency convolution of couplings and single loops (Peierls 
and Cooper) which decay as uj~‘^ [Eqs. (14) and 0]’ 
these high frequency singularities have in the end a lim¬ 
ited impact on susceptibilities and then correlations; the 
corrections remaining in most cases at a quantitative 
level. In special regions of the phase diagram, however, 
the consequence can be more substantial and can act on 
the nature of the ground state as will be discussed below. 

If we now move on the opposite side of the bisect¬ 
ing line, on the Hubbard side of the phase diagram at 
U = At, the strongest singularity is found for the SDW 
susceptibility at £c — 2.7, which comes with subdomi¬ 
nant singular BOW correlations [Fig. [^(b)]. These are 
the characteristics of the repulsive Hubbard model at 
half-filling which presents a charge or Mott gap A at 
£c, in accord with singular repulsive ^ 2 ( 0 , 0 , 0 ) and umk- 
lapp ^ 3 ^, 0 , 0 ) scatterings at zero frequency as meant 
in Fig. (b). As for the backseattering, the same Fig¬ 
ure shows that it becomes repulsive and large close to £c 
over the whole frequency range, with only a minimum for 
^ 1 ( 0 , 0,0) 5 ^ 2 , 3 ( 0 ,0, 0)]. Nevertheless, relevant repul¬ 

sive backscattering introduced by the dynamics leaves the 
spin sector gapless as shown by the singular behavior of 
the SDW response in Fig. [^(b). One also notes that the 
strong repulsive coupling that characterizes the backscat¬ 
tering amplitude gi{ooi,uj 2 ^oJs) at finite frequency has a 
limited influence on correlations, reducing for example 
the amplitude of BOW correlations compared to the sit¬ 
uation of the pure Hubbard model. 

Regarding the amplitude of the gap, the finite attrac¬ 
tive contribution coming from phonons [Eq. <§] weak¬ 
ens the amplitude of A compared to the pure Hubbard 
resulfP^, Ap = which yields £p = 7r/2 < £c, at 

U = At. 


1. Intermbehaviorediate phase 

As a function of U, the ground state evolves from the 
fully gapped CDW to the M-SDW states discussed above. 
Within the CDW state the reduction of the gap A with 
U carries on up to a critical value where static umklapp, 
^ 3 ( 0 , 0 , 0 ), is no longer sin gula r, but goes to zero and 
becomes irrelevant (Fig. whereas (0,0,0) and 

^ 2 ( 0 , 0,0) remain attractive as in Fig. (a). The irrel¬ 
evance of the static gs is concomitant with the onset of 
an additional, but sub dominant singular susceptibility in 
the SS channel, as shown in Fig. (c). From the same 
figure, we also note that the BOW response is also en¬ 
hanced. The CDW-SS singular correlations persist over 
a finite U -interval before the system ultimately enters in 
the M-SDW phase described above [Fig. (b)]. A de¬ 
tailed scan in the (U, ^fph)-plane allows us to outline an 
entire intermediate (I) region of the phase diagram of 
Fig. i (b) where similar features for CDW and SS cor- 
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FIG. 3: (Color online) Typical contour plot of umklapp 

scattering amplitude ^ 3 (ct;i,Ct; 2 , Ct;i) in the (oji, Ct; 2 )-plane close 
to in the intermediate I-phase at c^o = t. 



U/t 


relations are found. Note that within the I region, we 
observe a reinforcement of the SS correlations as U/t de¬ 
creases; its singular behavior reaching the CDW one as 
U ^ 0 and cjq ^ which is essentially the situation of 
the attractive Hubbard model. 

The boundary of the I state closes at some point 
whose locus depends on the phonon frequency 
ujQ. The resulting phase diagrams shown in Fig. for 
different ujq stand particularly well the comparison with 
previous numerical analysis, notably those carried out 
by the QMCP^ and DMRCt^ methods. Comparing for 
instance the metallic I region obtained from the QMC 
simulations of Hardikar and Clay (see Fig. 8 of Ref.l^^ 
with the one deduced here by RG, a quite accurate match 
is found for an amazing range of interactions and phonon 
frequencies. 

As to the the origin of the I phase, it must be stressed 
at the outset that within the electron gas model, irrele¬ 
vant umklapp and strong attractive backward scattering, 
along with degenerate CDW and SC singular suscepti¬ 
bilities, are well known characteri stics of the attractive 
ID Hubbard model at half fillin^^^^^^^. For the HH 
model, this limit is clearly realized at [/ = 0, along the 
^fph axis as cjo ^ 00 (see note in RefP). As ujQ is in¬ 
creased beyond t. Fig. shows indeed an ever-growing 
I region which includes the (0,^ph) axis. This suggests 
that the whole r egion is governed by a Luther-Emery 
type of fixed pointP^®, characterized by effective attrac¬ 
tive couplings and a gap in the spin sector, despite finite 
retardation and repulsive Both split the degeneracy 


FIG. 4: (Golor online) One-loop RG phase diagram of the 
ID Holstein-Hubbard model at different molecular phonon 
frequencies: ujq = 0.5t (a), ujq = t (b) and ujq = 5t (c). 


between CDW and SS correlations reinforcing the former 
with respect to the latter. According to the combination 
of couplings (18) for CDW susceptibility, singular attrac¬ 
tive umklapp scattering at finite frequency does have an 
impact, thought limited, in favoring an increase of CDW 
correlations and lack of degeneracy between CDW and 
SS in the I region [Fig. (c)]. 


The above results are compatible with those previously 
found from fRG by Tam et al^ for selected points of the 
phase diagram well outside and inside the supposed I re¬ 
gion. A different view was held, however, as to the prop¬ 
erties of charge degrees of freedom in this specific part of 
the phase diagram. The presence of singular and nega¬ 
tive umklapp scattering at finite frequency found in the 
I region (Fig. [^, was interpreted as the driving force of 
the CDW state, whose charge sector was then considered 
still gapped and insulating. We consider that a spatially 
uniform charge gapped state, if it exists, must be mani¬ 
fest in the equilibrium -static- properties rather than in 
the dynamics. Otherwise the Mott-insulating state thus 
obtained would be hard to reconcile with a finite charge 
compressibility,!!^ and singular superconducting correla¬ 
tions. 
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FIG. 5: (Color online) The power-law behavior 

Xbow ~ (straight lines) of the BOW susceptibil¬ 

ity on the LL line of the phase diagram at different U for the 
HH model in (a) [Fig. (b)], and for the PH model in (b) 
[Fig.|9] (c)]. In both cases, ojo/t = 1. 


2. Structure of the phase boundary 



We close this subsection by a detailed examination of 
the boundary to the M-SDW state in the phase diagram. 
As one moves along the bisecting line, ^ph = the RG 
calculations reveal the existence of a gapless transition 
line between the M-SDW and I phases (Fig. [^. At small 
U along this line, no singularity is encountered at finite £ 
in the susceptibilities. These rather exhibit a power law 
singularity of the form x ^ [Eo{£)]~^ at large £ as shown 
in Fig. (a). The exponent 7 is non universal, thought 
very small and positive, and is the largest for the BOW 
susceptibility. Within one-loop approximation, such a 
power law is characteristic of a Luttinger liquid (LL) with 
effective very weak repulsive interactions at low energy. 
This is confirmed by the contour plots of Fig. for the 
couplings in the (co’i,co’ 2 ) plane, where at large £ all the 
couplings are vanishingly small at the origin and remain 
weak at finite frequency. The LL line terminates on the 
bisecting line at the finite value t/* = which increases 
with cJo (Fig. [^. It has been checked that in the limit 
of large cjo, the LL and U = ^ph lines merge over the 
entire range of couplings with 7 -^ 0 , consistently with 
vanishing initial couplings in (|^ . It is worth noting that 
the existence of a LL metallic phase with a similar ujq 
dependence, albeit on a larger area of the phase diagram 
has been already pointed out by Fehske et using 
the DMRG technique, thought constrained by possible 
finite-size effects. 

Another feature emerges at finite ujq where the ending 
point (L'*,^ph) ill Fig. 0 marks the beginning of a dif¬ 
ferent boundary with the M-SDW phase. For U > EE 
a singularity appears in the susceptibilities at finite £c^ 
the strongest being for BOW, whereas SDW correlations 
become subdominant [Fig. (d)]; no enhancement in 


FIG. 6: (Golor online) Gontour plots of the scattering ampli¬ 
tudes in the (oji, 0 ^ 2 )-plane at large £ and uq — t 

for the metallic Luttinger liquid (LL) line of the HH model. 


CDW correlations is found. Actually, as one hits the 
boundary from below the relative importance of SDW 
and BOW correlations is inverted compared to the M- 
SDW phase [see Figs. (b), (d)]. The frequency pro¬ 
file of coupling constants in this BOW phase along the 
boundary is of interest. For small ^ph/^ and finite cjo, as 
shown in Fig. [!l( a), both 5^2(0, 0, 0) and ^3(0, 0, 0) scale to 
strong repulsive coupling as ^ ^ which is indicative of 
a charge gap. However, remains relatively 

small, although attractive at zero frequency and its close 
vicinity. It is from these frequency effects, a consequence 
of retardation, that comes the origin of dominant BOW 
correlations. According to the combination of couplings 
for the BOW susceptibility in Eq. (18), a change of sign in 
the backscattering in the low frequency range is sufficient 
to reinforce BOW correlations against SDW [Eq. (20)]. 

As g^Y^/t and U become larger at finite ujq along the 
boundary, 5 'i(co’i,Co’ 2 ,cji) develops much stronger attrac¬ 
tion at the origin and beyond, suggesting the presence 
of a spin gap at £c. This is reflected by the absence of 
singular SDW correlations in Eig. (d) (insert). While 
^ 2 ( 0 , 0,0) becomes small as shown in Eig. 0 (b), 53 re- 
mains positive at the origin and its neighbourhood, which 
indicates that at sizable ^ph the BOW order is now 
gapped in the spin sector and prevails over CDW that 
also becomes singular on the boundary. By cranking up 
further g^h/t at finite ujq along the boundary, namely be¬ 
yond the limitations of the RG, it is likely that BOW 























FIG. 8: (Color online) Selected susceptibilities versus £ in 
the P-BOW (a) and M-SDW (b) parts of the phase diagram 
of the Peierls-Hubbard model of Fig. [^(a) at ud = t. 


FIG. 7: (Color online) Contour plots of the scattering am¬ 
plitudes gi{(j 0 i,(j 02 ,(^i) in the (ct;i,Ct; 2 )-plane for the HH model 
in the BOW regime of the boundary with the intermediate 
phase at small U (a) and with the CDW phase at larger U 
(b). 


would be suppressed and the system would evolve to¬ 
ward a CDW ground state. The BOW phase emerging 
in this part of the CDW-M-SDW boundary, as a conse¬ 
quence of retardation at finite ujq, is reminiscent of the 
repulsive ID extended Hubbard lattice model where the 
BOW phase is known to enfold the U = 2V lin^^. This 
separatrix is known to separate the same CDW and M- 
SDW ground states in the continuum approximatiorP^. 
The ID extended Hubbard model being defined on a lat¬ 
tice, however, it is the momentum dependence of cou¬ 
plings, though irrelevant in the RG sense, that bre ak the 
CDW—M-SDW degeneracy in favor of a BOW stat d^^ l ^^ 


B. Peierls-Hubbard model 


1. V=0 case 


We now examine the PH model at finite repulsive U. 
As mentioned in Sec. [TTj the coupling of electronic bond 
transfer to acoustic phonons introduces particular ini¬ 
tial conditions for the phonon-mediated contribution to 
couplings of the PH model in Only the phonon 

induced backscattering is attractive, whereas the umk- 
lapp term, like 7/, is repulsive. According to (18), such a 
combination clearly favours the occurrence of the Peierls 
BOW (P-BOW) phase against CDW whose enhancement 
is totally absent across the (7/,^ph)—plane. As shown in 
Fig.i either a BOW or a SDW singularity is found at 
finite ^ph, as a function of U ; their tracking at H = 0 
yields the phase diagram of Fig. (top) at different De¬ 


bye frequencies. 

On the P-BOW side, the gap amplitude A = EFe~^^ 
extracted at £c [Fig. M (a)], is always larger than the 
U = 0 adiabatic mean-field value, Aq = Fi¬ 

nite U enhances BOW correlations and then s treng thens 
their coupling to 2k f phonons in Eqs. (The 
enhancement of the Peierls order parameter by repulsive 

is a well-known result of early numerical simulations 
on that model^. Regarding the scattering amplitudes, 
singularity of ( 0 , 0 , 0 ) is found in the attractive sec¬ 
tor, while the strong-coupling peaks of both ^2 ( 0 , 0 , 0 ) 
and umklapp ^ 3 ( 0 , 0 , 0 ) are for positive values (Fig. [Tq|). 
These are indicative of a fully gapped phase, as expected 
for the commensurate Peierls order at half-filling. The 
frequency dependence of the gi{uJi, 002 ^^ 1 ) displays simi¬ 
lar features over the whole P-BOW region. However, the 
influence of strong coupling at finite frequency on corre¬ 
lations remains weak. 

The transition line between P-BOW and M-SDW 
states in the phase diagram is also of interest. Owing to 
the reinforcement of electronic BOW correlations by 7/, 
the transition to the M-SDW line occurs below the bisect¬ 
ing line ^ph = 7/, in contrast to what has been obtained 
for the HH model (Fig.[^. As a result, moving toward the 
adiabatic limit for small ujd the P-BOW phase grows in 
importance to the detriment of M-SDW. In Fig. (bot¬ 
tom), we have plotted the phonon coupling parameter 
a[= {gphOJD/SCV^'^] versus U for the transition linM 
The trace fits fairly well the power-law ac ^ , with 

the exponent 77:^0.31atR = 0. Similar algebraic varia¬ 
tion has been rep orted by QMC simulations at finite 
(see Sec. HIB 2 | ). For most ^fph values, the line ac{U) 
shows a direct transition between the P-BOW and M- 
SDW phases [Fig. (top)], except for very small ^ph, 
where there is a finite but minute U interval in which LL 
metallic conditions prevail (see insets of Fig. |^. As in the 
HH model, a power-law dependences of the BOW suscep¬ 
tibility is found, as shown in Fig. [^(b). Otherwise, due 
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FIG. 9: (Color online) Phase diagram of the PH model for 
(top) ujolt — 0.1,1, and 5 at H = 0; (middle) — 1 

and y > 0 (inserts: zoom of the Luttinger liquid line at 
very small couplings); (bottom) comparison with a power law 
dependence of the boundary line between P-BOW and M- 
SDW (dashed lines) for different values of V (bottom). The 
full circles are the QMC results of Ref.^. 


to the fact that for the PH model both the Coulomb and 
phonon-induced terms contribute positively to umklapp 
scattering [Eqs. © and ([Tq|)], the charge sector is gapped 
and the conditions for the emergence of an intermediate 
metallic phase are not satisfied. 

As to the M-SDW phase at sufficiently strong [/, their 
characteristics are similar to the one found for the HH 
model. This is confirmed by the frequency profiles of 
scattering amplitudes of Fig. (b) close to when 
compared to those of Fig.|^(b). The backscattering am¬ 
plitude develop sizable repulsive values with a minimum 
for ^i(0, 0, 0), subordinated to the singular positive peaks 
^ 2 ( 0 , 0,0) and ^ 3 ( 0 , 0, 0) at positive values. The Mott gap 
is associated with the strongest singularity in xsdw at ic 
closely followed by Xbow, as shown in Fig. |(b). 


2. Finite V 

We close the section by examining the influence of 
small nearest-neighbor repulsion V on the structure of 



FIG. 10: (Golor online) Gontour plot of the scattering ampli¬ 
tudes gi(cJi, UJ 2 , cui) in the (oji, 0 ^ 2 )-plane for the Peierls Hub¬ 
bard phase diagram of Fig. (top) close to £c and at cud = t 
: (a), the P-BOW phase and (b), the M-SDW phase. 


the phase diagram of the PH model. In Fig. (middle), 
the critical line ac{U) is plotted for different V (< U/2). 
We first observe that the P-BOW region increases in size 
with y, as described by the fit to the expression ac = aU'^ 
showing a drop of the coefficient a and a power-law ex¬ 
ponent T] which depends on the value of V. The rein¬ 
forcement of BOW correlations and in turn of the Peierls 
state follows from the decrease (increase) of gi {g 2 ) by 
y in Eq. (10), which compensates the drop in umklapp 


and then according to ( [l^ makes the BOW correlations 
and its coupling to phonons larger. The frequency depen¬ 
dence of the ^i(co’i, 0 ^ 2 , cJi) projected in the (c^i, Co’ 2 )-plane 
present essentially the same characteristics as those ob¬ 
tained in the Hubbard case for y = 0 (Eig. [1Q| ). 

The results obtained at y = U/A can be com¬ 
pared to QM C simulations performed in the same model 
conditionsP^ As shown in Eig. [^(bottom), the power law 
found by QMC for ac is fairly well reproduced by the RG 
method over the whole interval of acceptable couplings 
for a weak coupling scheme. 


IV. SUMMARY AND CONCLUDING 
REMARKS 

Summarizing, we have studied the phase diagrams 
of ID Holstein-Hubbard and Peierls-Hubbard models at 
half filling by a weak-coupling RG method. The {U 
phase diagrams that are mapped out from the susceptibil¬ 
ities and frequency-dependent scattering amplitudes at 
the one-loop level adhere for the most part to the results 
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of numerical simulations, in particular those of quantum 
Monte Carlo methods for which the most detailed stud¬ 
ies have been carried out and the best comparison can 
be established in weak coupling. 

In the case of the Holstein-Hubbard model, the RG re¬ 
sults pinpoint a precise delimitation of an intermediate 
region separating the commensurate CDW and a Mott- 
insulating SDW phases in the phase diagram. The sin¬ 
gular static superconducting correlations found through¬ 
out this region, though always subordinate to CDW and 
tied to irrelevant umklapp scattering at zero frequency, 
are congruent with the metallic character of this region. 
The present results also revealed an internal structure 
of the phase boundary with the M-SDW phase with 
the presence of distinct phases, unanticipated and un¬ 
explored from the viewpoint of previous works on this 
model. A very weakly correlated Luttinger liquid fol¬ 
lowed by a BOW phase characterized by either a charge 
or spin gap has been found depending on the strength of 
the Coulomb term. The BOW phase replaces the Mott- 
SDW state as a result of retardation, whose dynamics 
generates an attraction in the backscattering amplitude 
alone, in analogy with momentum-dependent effects for 
the appearance of the BOW phase enfolding the U = 2V 


line of the ID extended Hubbard lattice model. 

A similar analysis carried out on the Peierls-Hubbard 
model with its peculiar momentum- and frequency- 
dependent scattering amplitudes rather reveals a Mott 
SDW state competing with the Peierls-BOW order whose 
prominence in the phase diagram grows as a function of 
retardation, and intrasite and inter site Coulomb terms 
U and V. The nearly structureless boundary between 
the two ground states in the (t/, ^ph)—plane of the phase 
diagram follows a power-law profile compatible with nu¬ 
merical simulations. 
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